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Abstract 

We quantify the impact of gauge anomalies at the Large Hadron Collider by studying the invariant 
mass distributions in Drell-Yan and in double prompt photon, using an extension of the Standard Model 
characterized by an additional anomalous U{\) derived from intersecting branes. The approach is rather 
general and applies to any anomalous abelian gauge current. Anomalies are cancelled using either the 
Wess-Zumino mechanism with suitable Peccei-Quinn-like interactions and a Stiickelberg axion, or by the 
Green-Schwarz mechanism. We compare predictions for the corresponding extra Z prime to anomaly-free 
realizations such as those involving U{1)b-l- We identify the leading anomalous corrections to both 
channels, which contribute at higher orders, and compare them against the next-to-next-to-leading order 
(NNLO) QCD background. Anomalous effects in these inclusive observables are found to be very small, far 
below the percent level and below the size of the typical QCD corrections quantified by NNLO X-factors. 
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1 Introduction 



The study of anomalous gauge interactions at the LHC and at future hnear cohiders is for sure a difficult 
topic, but also an open possibility that deserves close theoretical and experimental attention. Hopefully, 
these studies will be able to establish if an additional anomalous extra Z' is present in the spectrum, 
introduced by an abelian extension of the gauge structure of the Standard Model (SM), assuming that extra 
neutral currents will be found in the next several years of running of the LHC [T] . The interactions that we 
discuss are characterized by genuine anomalous vertices in which gauge anomalies cancel in some non trivial 
way, not by a suitable (anomaly-free) charge assignment of the chiral fermion spectrum for each generation. 

The phenomenological investigation of this topic is rather new, while various mechanisms of cancellation 
of the gauge anomalies involving axions have been around for quite some time. Global anomalies, for 
instance, introduced for the solution of the strong CP problem, such as the Peccei-Quinn (PQ) solution 
[SlEHHEKlllTllHllQ] (reviewed in [10]) require an axion, while local anomalies, cancelled by a Wess-Zumino 
counterterm, allow an axion-like particle in the spectrum, whose mass and gauge coupling - differently from 
PQ axions - are independent. Similar constructions hold also in the supersymmetric case and a generalization 
of the WZ mechanism is the Green-Schwarz mechanism (GS) of string theory. The two mechanisms are 
related but not identical [11], the first of them being characterized by a unitarity bound. Details on the 
relation between the two at the level of effective field theory can be found in [IH I12j . 

Intersecting brane models, in which several anomalous C/(l)'s and Stiickelberg mass terms are present, 
may offer a realization of these constructions [131 [T^ [T5\ 116] . which can also be investigated in a bottom- 
up approach by using effective lagrangeans built out of the requirements of gauge invariance of the 1-loop 
effective action |17] . In our analysis we will consider the simplest extension of these anomalous abelian gauge 
factors, which involves a single anomalous U{1), denoted as U{1)b- The corresponding anomalous gauge 
boson (B) gets its mass via a combination of the Higgs and of Stiickelberg mechanisms. Axions play a key 
role in the cancellation of the anomalies in these theories although they may appear in other constructions 
as well, due to the decoupling of a chiral fermion in anomaly-free theories [18] . 

The presence of an anomalous U{1) in effective models derived from string theory is quite common, 
although in all the previous literature before [T7] and [18|. \T9\ [20] the phenomenological relevance of the 
anomalous U{1) had not been worked out in any detail. In particular, the dynamics of the anomalous extra 
gauge interaction had been neglected, by invoking a decoupling of the anomalous sector on the assumption 
of a large mass of the extra gauge boson. In [17j it was shown that only one physical axion appears in 
the spectrum of these models, independently of the number of abelian factors, which is the most important 
feature of these realizations. In our case, the axion can be massless or massive, depending on the structure 
of the scalar potential. Recent developements in the study of these models include their supersymmetric 
extensions [21] and their derivations as symplectic forms of supergravity [221 123] . Other interesting variants 
include the Stiickelberg extensions considered in [24:\ [25l [26] which depart significantly from the Minimal 
Low Scale Orientifold Model (MLSOM) introduced in [T7| and discussed below. Specifically these models 
are also characterized by the presence of two mechanisms of symmetry breaking (Higgs and Stiickelberg) 
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but do not share the anomalous structure. As such they do not describe the anomalous C/(l)'s of these 
special vacua of string theory. 

Axion-like particles, beside being a natural candidate for dark matter, may play a role in explaining 
some puzzling results regarding the propagation of high energy gamma rays [27\ |28] due to the oscillations 
of photons into axions in the presence of intergalactic magnetic field. In general, the presence of independent 
mass/coupling relations for these particles allows to evade most of the experimental bounds coming from 
CAST and other experiments on the detection of PQ axions, characterized by a suppression of both mass 
and gauge couplings of this particle by the same large scale (10^° — 10^^ GeV), (see [291 EQ]). While a 
phenomenological study of Stiickelberg axions is underway in a related work, here we focus our attention on 
the gauge sector, quantifying the rates for the detection of anomalous neutral currents at the LHC in some 
specific and very important channels. 

• Drell-Yan 

Being leptoproduction the best way to search for extra neutral interactions, it is then obvious that the study 
of the anomalous vertices and of possible anomalous extra Z' should seriously consider the investigation 
of this process. We describe the modifications induced on Drell-Yan computed in the Standard Model 
(SM) starting from the description of some of the properties of the new anomalous vertices and of the 
corresponding 1-loop counterterms, before moving to the analysis of the corrections. These appear - both in 
the WZ and GS cases in the relevant partonic channels at NNLO in the strong coupling constant (0(a^)). 
We perform several comparisons between anomalous and non-anomalous extra Z' models and quantify the 
differences with high accuracy. 

• Direct Photons (Di-photon, DP) 

Double prompt (direct) photons offer an interesting signal which is deprived of the fragmentation contribu- 
tions especially at large values of their invariant mass Q, due to the steep falling of the photon fragmentation 
functions. In addition, photon isolation may provide an additional help in selecting those events coming 
from channels in which the contribution of the anomaly is more sizeable, such as gluon fusion. Also in this 
case we perform a detailed investigation of this sector. For direct photons, the anomaly appears in gluon 
fusion -at parton level - in a class of amplitudes which are characterized by two-triangles graphs - or BIM 
amplitudes - using the definitions of [11] . 

In both cases the quantification of the background needs extreme care, due to the small signal, and the 
investigation of the renormalization/factorization scale dependence of the predictions is of outmost impor- 
tance. In particular, we consider all the sources of scale-dependence in the analysis, including those coming 
from the evolution of the parton densities (Pdf's) which are just by themeselves enough to overshadow 
the anomalous corrections. For this reason we have used the program Candia in the evolution of the 
Pdf's, which has been documented in [31j. The implementations of DY and DP are part of two programs 
CANDiA£)y and C ANDiA Axion for the study of the QCD background with the modifications induced by the 
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anomalous signal. The QCD background in DP is computed using Diphox |32| and Gamma2MC based 
on Ref. [33]. The NLO corrections to DP before the implementation of Diphox have been computed by 
Gordon and one of the authors back in 1995 [M] and implemented in a Monte Carlo based on the phase 
space slicing method. 

In the numerical analysis that we present we have included all the contributions coming from the two 
mechanisms as separate cases, corrections that are implemented in DY and DP. We will start analyzing the 
contributions to these processes in more detail in the next sections, discussing the specific features of the 
anomalous contributions and of the corresponding counterterms at a phenomenological level. To make our 
treatment self-contained we have summarized some of the properties of effective vertices in these models, 
relevant for our analysis. 

Our work is organized as follows. After a brief description of the anomalous interactions and of the 
counterterms that appear either at lagrangean level (WZ case) or at the level of the trilinear gauge vertex 
(GS case), we discuss the main properties of these vertices and we address the structure of the corrections in 
DY and in DP. Our study is mainly focused on the invariant mass distributions in the two cases. The need 
for performing these types of analysis in parallel will be explained below, and there is the hope that it may 
be extended to other processes and observables in the future, such as rapidity distributions and rapidity 
correlations [35]. We present high precision estimates of the QCD background at NNLO, which is the order 
where, in these processes, the anomalous corrections start to appear. Other analysis, of course, are also 
possible, such as those involving 4-fermion decays in trilinear gauge interactions which could, in principle, be 
sensitive to Chern-Simons terms [171I36| if at least two anomalous C/(l)'s are present in the spectrum. These 
additional interactions are allowed [37j whenever the distribution of the partial anomalies on a diagram is 
not fixed by symmetry requirements. A complete description of these vertices has been carried out in |37j . 
useful for direct phenomenological studies. This search, however, is expected to be experimentally also very 
difficult. 

As we are going to show, the search for effects due to anomalous U{1) at the LHC in pp collisions cannot 
avoid an analysis of the QCD background at NNLO. DY and DP are the only two cases where this level 
of precision has been obtained in perturbation theory. As we are going to show, the anomalous effects at 
the LHC in these two key processes are tiny, since the invariant mass distributions are down by a factor of 
10^ - 10^ compared to the NNLO (QCD) background. The accuracy required at the LHC to identify these 
effects on these observables should be of a fraction of a percent (0.1% and below), which is beyond reach at 
a hadron collider due to the larger indetermination intrinsic in QCD factorization and the parton model. 

2 Anomaly-free versus an anomalous extra Z' in Drell-Yan 

As we have already mentioned, the best mode to search for extra Z' at the LHC is in the production of 
a lepton pair via the Drell-Yan mechanism (qq annihilation) mediated by neutral currents. The final state 
is easily tagged and resonant due to the s-channel exchange of the extra gauge boson. In particular, a 



4 



new heavier gauge boson modifies the invariant mass distribution also on the Z peak due to the small 
modifications induced on the couplings and to the Z — Z' interference. In the anomalous model that we 
have investigated, though based on a specific charge assignment, we find larger rates for these distributions 
both on the peak of the Z and of the Z' compared to the other models investigated, if the extra resonance 
is around 1 TeV. This correlation is expected to drop as the mass of the extra Z' increases. In our case, as 
we will specify below, the mass of the extra resonance is given by the Stiickelberg (Mi) mass, which appears 
also (as a suppression scale) in the interaction of the physical axion to the gluons and is essentially a free 
parameter. 

In DY, the investigation of the NNLO hard scatterings goes back to j[38j , with a complete computation of 
the invariant mass distributions, made before that the NNLO corrections to the DGLAP evolution had been 
fully completed. In our analysis we will compare three anomaly-free models against a model of intersecting 
brane with a single anomalous f/(l). The anomaly- free charge assignments come from a gauged B — L 
abelian symmetry, a "g -|- u" model -both described in ^9] - and the free fermionic model analyzed in |40j . 
We start by summarizing our definitions and conventions. 

In the anomaly-free case we address abelian extensions of the gauge structure of the form SU{3) x 
SU{2) X U{1)y X [7(1)2, with a covariant derivative in the Wj^, By, Bz (interaction) basis defined as 



where we denote with g2,gY,9z the couplings of SU{2), U{1)y and U{1) 
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where e is a perturbative parameter which is around 10 ^ for the models analyzed, introduced in |39j and 
lOl. It is defined as 
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In this class of models we have two Higgs doublet Hi and H2, whose vevs are vhi and VH2 and an extra 
SU{2)\Y singlet whose vev is v^. The extra U{l)z charges of the Higgs doublet are respectively and 
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zh2, while for the singlet this is denoted as z^. Taking the value of of the order of the electroweak scale 
(~ 246 GeV), we fix vhi with tan/? = vh2/'vHi, and we still have one free parameter, w^, which enters in the 
calculation of the mass of the extra Z'. Then it is obvious that we can take the mass Mz' and the coupling 
constant as free parameters. We choose tan /3 « 40 in order to reproduce the mass of the Z boson at 
91.187 GeV, choice that is performed, for consistency, also in the anomalous model. In this last case the 
Higgs sector is characterized only by 2 Higgs doublets, with the vev of the extra singlet being replaced by 
the Stiickelberg mass. We define §2 sin 6w = Qy cos 9w = e and construct the charge eigenstates and 
the corresponding generators as usual 
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which relates the interaction and the mass eigenstates. Substituting these expression in the covariant 
derivative we obtain 
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where we have neglected all the O(e^) terms. Sending gz ^ and e — > we obtain the SM expression. 
The vector and the axial couplings of the Z and Z' to the fermions are expressed equivalently in terms of 
the left - (zl) and right - (zr) U{1)z chiral charges and hypercharges (Y^j, Yl) of the models that we have 
implemented. These can be found in [40J for the free fermionic case and in [3U] for the remaining models 
with a V-A structure given by 
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where j is an index which represents the quark or the lepton and we have set sin Q-^^ = , cos 9w = Cw for 
brevity. 



2.1 An anomalous extra Z' 

In presence of anomalous interactions we can use the same formahsm developed so far for anomaly-free 
models with some appropriate changes. Since the effective lagrangean of the class of the anomalous models 
that we are investigating includes both a Stiickelberg and a two-Higgs doublet sector, the masses of the 
neutral gauge bosons are provided by a combination of these two mechanisms. In this case we take as 
free parameters the Stiieckelberg mass Mi and the anomalous coupling constant qb, with tan/3 as in the 
remaining anomaly-free models. As we have already stressed, the analysis does not depend significantly on 
the choice of this parameter. The value of the Stiickelberg mass Mi is loosely constrained by the D-brane 
model in terms of suitable wrappings (n) of the 4-branes which define the charge embedding |4H ll4j reported 
in Tabs. EEElandgl 

The mass-matrix in the neutral gauge sector is given by 
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where B is the Stiickelberg field and the mass matrix is defined as 
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with 



Nbb = {Qu'^Vu + qPvI) 9I Xb = {q^vl + q^vj) 9b- (10) 

Here and denote the vevs of the two Higgs fields while and q^ are the Higgs charges 

under the extra anomalous \J(X)b- We have also defined v = \j^^^t^r^ and g = y'^si + fi'y • massless 
eigenvalue of the mass matrix is associated to the photon j4^, while the two non-zero mass eigenvalues 
denote the masses of the Z and of the Z' vector bosons. These are given by 
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- 2 "Mf 4 +M4 8 ^^BB-gv), 
Ml, = i (^2Mf + 5^^;' + A^BB + V (2M2 - g^v^ + iV^^)' + ig^xj,^ (12) 
^ 2 

The mass of the Z gauge boson gets corrected by terms of the order v'^/Mi, see Fig. [H converging to the 
SM value as Mi — > cxd, while the mass of the Z' gauge boson can grow large with Mi. The physical gauge 
fields can be obtained from the rotation matrix 
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which can be approximated at the first order as 
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which is the analogue of the matrix in Eq. ^ for the anomaly-free models, but here the role of the mixing 
parameter ei is taken by the expression 



XB_ 

A relation between the two expansion parameters can be easily obtained in an approximate way by a direct 
comparison. For simplicity we take all the charges to be 0(1) in all the models obtaining 
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Ml-Ml^gy^ 
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giving 



which is the analogue of Eq. p5]) , having identified the Stiickelberg mass with the vev of the extra singlet 
Higgs, Ml ~ gzV^. This is natural since the Stiickelberg mechanism can be thought of as the low energy 
remnant of an extra Higgs whose radial fluctuations have been frozen and with the imaginary phase surviving 
at low energy as a CP-odd scalar 
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Figure 1: Anomalous A4z' as a function of tlie coupling gs- 
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Table 1: Parameters for a Class A model with a D6-brane. 

Concerning the charge assignments, the corresponding model is obtained form the intersection of 4 
branes (a, b, c, d) with generators {qa, qt, qc, Qd) which are rotated to the hypercharge basis, with an anomaly 
free hypercharge. The U{l)a and ^7(l)d symmetries are proportional to the baryon number and the lepton 
number respectively. The U{l)c symmetry can be considered as the third component of the right-handed 
weak isospin; the U{l)b is a PQ-like symmetry. A discussion of the construction can be found in |14j and 
|41j . The identification of the generators involve the solution of some constraint equations, solutions which 
for a compactification are parametrized by a phase e = ±1; the Neveu-Schwarz background on the first 
two tori /3j = 1 — 6j = 1, 1/2, four integers na2,nbi,nci,nii2, which are the wrapping numbers of the branes 
around the extra (toroidal) manifolds of the compactification, and a parameter p = 1, 1/3, with an additional 
constraint in order to obtain the correct massless hypercharge. One of the choice for these parameters is 
reported in Table [TJ 

3 The GS and WZ vertices and gluon fusion 

As we have mentioned above in the previous sections, the two available mechanisms that enforce at the level 
of the effective lagrangean the cancellation of the anomalies involve either PQ-like (axion-like) interactions 
- in the WZ case - or the subtraction of the anomaly pole (for the GS case). In a related analysis [12] we 
have presented some of the main features of the two mechanisms taking as an example an axial (anomalous) 
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Table 2: Higgs charges in the Madrid modeL 
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Table 3: SM spectrum charges in the D-brane basis for the Madrid model. 

version of QED to illustrate the cancellation of the anomaly in the two cases. In the GS case, the anomaly 
of a given diagram is removed by subtracting the longitudinal pole of the triangle amplitude in the chiral 
limit. We have stressed in [S^ that the counterterm (the pole subtraction) amounts to the removal of one of 
the invariant amplitudes of the anomaly vertex (the longitudinal component) and corresponds to a vertex 
re-definition. 

The procedure is exemplified in Fig. [2] where we show the triangle anomaly and the pole counterterm 
which is subtracted from the first amplitude. The combination of the two contributions defines the GS 
vertex, which is made of purely transverse components in the chiral limit [T^] and satisfies an ordinary Ward 
identity. Notice that the vertex does not require an axion as an asymptotic state in the related S-matrix; 
for a non-zero fermion mass in the triangle diagram, the vertex satisfies a broken Ward identity. We now 
proceed and summarize some of these properties, working in the chiral limit. 

Processes such as gg 77, mediated by an anomalous gauge boson Z', can be expressed in a simplified 
form in which only the longitudinal component of the anomaly appears. We therefore set kf = k2 = and 
ruf = 0, which are the correct kinematical conditions to obtain the anomaly pole, necessary for a parton 
model (factorized) description of the cross section in a pp collision at the LHC, where the initial state of the 
partonic hard-scatterings are on-shell. 

We start from the Rosenberg form of the AVV amplitude, which is given by 

+A4k!^£[ki,k2,i^,X] + A5k'(e[ki,k2,fJ.,X] + A(ik^e[ki,k2, n, X] , (18) 
and imposing the Ward identities to bring all the anomaly on the axial-vector vertex, we obtain the usual 
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Table 4: Fermion spectrum charges in the y-basis for the Madrid model [IT] . 




Figure 2: A gauge invariant GS vertex of the AW type, composed of an AW triangle and a single counterterm of 
the Dolgov-Zakharov form. 



conditions 

Ai = klAi + ki ■ k2A3 
A2 = kfA^ + ki ■ feavle 
^3(^:1,^2) = -^6(^1,^2) 

^4(A;i,A;2) = -^5(^1,^2), (19) 

where the invariant amplitudes A^, . . . , Aq are free from kinematical singularities for off-shell external lines. 
We set k'^ = {ki + k2)^ = s. As we have mentioned, in the parton model we take the initial gluons to 
be on-shell, while the hadronic cross section is obtained by convoluting the hard scattering given above 
(corrected by a color factor) with the Pdf's. The amplitude simplifies drastically in this case and takes the 
form 

T^'^^ = AQk^e[ki,k2,y, ^A + (^4 + ^e) (A:^e[A;i, fca, ;U, A] - k^e[ki,k2, v,\]), (20) 

in which the second piece drops off for physical on-shell photon/gluon lines, leaving only a single invariant 
amplitude to contribute to the final result 

T^'^^ = 4{s){ki + k2)^e [ku k2, u, ^l] (21) 

where 

1 / 77i^ p -|- 1 \ / TnP' 
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Figure 4: Ward identity in the WZ case in the chiral limit. 



The anomaly pole is given by the first term of Eq. (j22p 

= -l-{ki + k2)^e[h,k2,iy,fi] . (23) 
The logarithmic functions in the expression above are continued in the following way in the various region 
< s < Amj : 



s > Amj > : 



Notice that the surviving amplitude Aq multiplies a longitudinal momentum exchange and, as discussed 
in the literature on the chiral anomaly in QCD 142^143] . is characterized by a massless pole in s, which is the 
anomaly pole, as one can clearly conclude from Eq. (j22p . This equation shows also how chiral symmetry 
breaking effects appear in this amplitude at this special kinematical point by the nif terms. 

The subtraction of the anomaly pole is shown in Fig.[2]and is represented by diagram c). The combination 
of diagrams b) and c) defines the GS vertex of the theory [12], with diagram c) described by Eq. ([23|) (— Tc). 
It is easily verified that in the massless fermion limit and for on-shell gluon lines, the GS vertex is trivially 
vanishing by construction. In general, for any asymmetric configuration of the external lines in the vertex. 
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Figure 5: Generalized Ward identity in the WZ case. 



even in the massless limit, the vertex has non-zero transverse components |441 05]. The expression is well 
known |461 145j in the chiral limit and has been shown to satisfy the Adler-Bardeen theorem |45j . 

For a non- vanishing mj the GS vertex, for generic virtualities, can be defined to be the general AVV 
vertex, for instance extracted from [47] or, in the longitudinal/transverse formulation, by the amplitudes 
given in [35], with the subtraction of the anomaly pole, as given in [12j. We will refer to the anomaly 
(subtraction) counterterm of diagram b) as to the Dolgov-Zakharov [33] (DZ) counterterm. The anomaly 
diagram reduces to its DZ form for two on-shell gauge lines (photons /gluons) and in this case the transverse 
components completely disappear. There are other cases in which, instead, the longitudinal components 
cancel. This occurs if, for instance, a conserved current is attached to the anomalous line, rendering the 
anomaly "harmless", as explained in [12]. 




The analogous interaction in the WZ case is shown in Fig. [3l where we have attached a fermion pair in the 
final state to better identify the contributions. In this case, beside the anomalous contribution of diagram 
a), the mechanism will require the exchange of a physical axion, shown in diagram b) and c). Diagram b) 
is the usual WZ counterterm (or generalized PQ interaction) while the third diagram is non-vanishing only 
in the presence of fermions of non-zero mass. This third contribution is numerically irrelevant and in DY 
is usually omitted. The WZ mechanism re-establish gauge invariance of the effective lagrangean but is not 
based on a vertex re-definition and, furthermore, involves an asymptotic axion state. As shown in jllj the 
presence of a unitarity bound in this mechanism is a signal of its limitation as an effective theory (see also 
the discussion in [12]). We have summarized in an appendix the discussion of this point in a simple case. 

3.1 Ward identities 

Both vertices satisfy ordinary Ward identities in the chiral limit and generalized Ward identities away from 
it. In the chiral limit, for instance, the WZ mechanism adds to the effective action of the anomalous theory 
an interaction of the Stiickelberg field {b) with the gluons (bG AG), shown in diagram b) of Fig. HI In this 
figure we have shown a diagrammatic realization of the Ward identity for this case. 

In WZ, being the cancellation based on a local field theory, the derivation of the generalized Ward 
identity can be formally obtained from the requirement of BRST invariance of the gauge-fixed effective 
action, as we have shown in [37]. This is illustrated in Fig. [5l in the case of an anomalous Z', where we 
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Figure 6: Generalized Ward identity in the GS case. 
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Figure 7: BIM amplitude with an anomalous Z' exchange. 



show the coupling of the goldstone - in the broken Higgs phase - to the gluons (diagram b)) and to the 
massive fermion (diagram c)) [37J. The normalization of the counterterm in b) can be chosen to remove 
the anomaly of diagram a) when a single fermion runs inside the anomaly loop. Alternatively, the same 
graphical representation holds if in the first and the last diagram we sum over the entire generation. In this 
case the counterterm is normalized to cancel the entire anomaly of the complete vertex. 

The analysis in the GS case is slightly more subtle. We show in Fig. [6] the generalized Ward identity 
satisfied by the vertex in the massive case. In the massless case only diagrams a) and b) survive, while 
the contribution of diagram c) comes from a direct computation. It is obtained by multiplying typical 
pseudoscalar interaction - such as the one shown in Fig. [5l diagram c) - by a massless pole. If we denote by 
T^'^ the diagram describing the decay of a pseudoscalar into two gluons, diag ram c) takes the form j Iz^T^y ^ 
with a factorized pole on the anomalous external line. We refer to [12] for a detailed discussion of these 
points. 

In our analysis we encounter a class of amplitudes (BIM ) [11] which are characterized by two anomaly 
vertices connected by an s-channel exchange of the anomalos gauge boson. These amplitudes grow quadrat- 
ically with the energy and are not eliminated by fine tuning Fig. [71 The true BIM amplitude is the one 
shown in diagram b) and appears in the gluon fusion sector in the WZ case. In the SM a similar graph 
contributes only if heavy fermions run in the loop. They are comparable in size to the anomalous BIM am- 
plitude. Obviously, this contribution would be identically vanishing if all the fermions of a given generation 
would be mass-degenerate. 

Diagram a) shown in the same figure, instead, is the GS version of the BIM amplitude and is identically 
zero in the chiral limit for on-shell gluon lines, as is the case in the parton model. For this reason the gluon 



14 



fusion sector disappears completely for DP (in the GS mechanism), since the BIM amplitude in this case is 
obtained by replacing diagram b) of this figure with diagram a) which is indeed vanishing. 

We can summarize the basic features of the anomalous sectors in anomaly free-models after QCD fac- 
torization, for generic virtualities of the external gauge lines, according to the following points: 

1) In the SM the residual contributions coming from anomalous diagrams, such as in the VVZ, VVZ' 
vertices, where ^ is a gauge field, are proportional to the mass of the heavy quarks in the anomaly loop. In 
the chiral limit, instead, both the anomaly pole contribution and the transverse component of the anomaly 
cancel by charge assignment. 

2) In the GS case, as we have just discussed, the anomaly pole is absent by definition, while the trans- 
verse contributions are allowed. This separation between longitudinal and transverse components is less 
transparent for a heavy fermion mass, which induce a longitudinal component, proportional to m^/s^ times 
a small logarithmic correction of the ratio of the same variables, away from the chiral limit. This longitu- 
dinal component, however, should not be confused with the anomaly pole and is not shifted or corrected 
perturbatively in any way. It can couple, for instance, to a tt (top) quark current because of a broken Ward 
identity and can be interpreted as a manifestation of the GS mechanism at the LHC, but can be easily 
overshadowed by SM contributions. This point will be re-addressed more formally below in Eq. (j32p . 

At some special kinematical points (two massless gauge lines, or three massless gauge lines of the same 
virtualities) where the anomalous vertex takes its DZ form, the GS vertex is identically vanishing in the 
massless case. In the presence of a heavy fermion the logarithmic correction shown in Eq. (j22p reappears. 

3) In the WZ case the anomaly pole is not cancelled. A second sector (the exchange of the axion) is 
needed to restore the gauge invariance of the effective action. In a hadronic collision the BIM amplitudes 
induce very small deviations from the SM behaviour after the convolution with the gluon density. They are 
absent in DY at NNLO. In DP they affect the invariant mass distributions - at large Q - of the photon pair, 
for a given center of mass energy of the two colliding protons. As such they are sensitive to large (Bjorken) 
x-values of the gluon Pdf's, region where the gluon density is rapidly decreasing. In particular, in DP their 
contribution becomes more sizeable via intereference with some box-like amplitudes {gg 77). In previous 
NLO study of this process [M] they had been included even though they exceed the NLO accuracy, being 
truly NNLO contributions. These amplitudes and vertices are the basic building blocks of our numerical 
analysis and are responsible for all the anomalous signal both in DY and in DP. We will try to quantify 
their impact in the invariant mass distributions in both cases. 

4 Invariant mass distributions in Drell-Yan 

Our NNLO analysis of the invariant mass distributions for lepton pair production, for the computation of the 
QCD sectors, is based on the hard scatterings of ^8], and the NNLO evolution of the parton distributions 
(Pdf's) has been obtained with Candia |31j . The anomalous corrections to the invariant mass distributions 
have been evaluated separately, since at NNLO they appear in DY in the interference with the lowest order 
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graph, and added to the standard QCD background. It is important to recall that lepton pair production 
at low Q via Drell-Yan is sensitive to the Pdf's at small-x values, while in the high mass region this process 
is essential in the search of additional neutral currents. In our analysis we have selected a mass of 1 TeV 
for the extra gauge boson and analyzed the signal and the background both on the peaks of the Z and of 
the of the new resonance. 

At hadron level the colour-averaged inclusive differential cross section for the reaction Hi + H2 — 
^1 + ^2 + is given by the expression 



raz{Q\Ml)Wz{r,Q^) r = ^, (25) 



where Z = Z,Z' is the point-like cross section and all the information from the hadronic initial state is 
contained in the Pdf's. The hadronic structure function Wz{t,Q'^) is given by a convolution product 
between the parton luminosities <^jj(x, and the Wilson coefficients Ajj(a;, Q^, /i^, 



■ ■ Jt X 



(26) 



where the luminosities are given by 



J X y \y J 

and the Wilson coefficients (hard scatterings) depend on both the factorization (/Ui?) and renormalization 
scales (^_r), formally expanded in the strong coupling as as 



A,,(x,Q2,;u2,) = f;<(/x|)A;;)(x,Q2,;u2„/,|). (28) 

n=0 



We will vary \ip and [Ir independently in order to determine the sensitivity of the prediciton on their 
variations and their optimal choice. 

The anomalous corrections to the hard scatterings computed in the SM will be discussed below. We just 
recall that the relevant point-like cross sections appearing in the factorization formula (j25p and which are 
part of our analysis include, beside the Z and the Z' resonance, also the contributions due to the photon 
and the 7 — Z, 7 — Z' interferences. For instance in the Z' case we have 



3Q4 



AMz' sin2 9w cos2 OwNc {Q^ - M|,)2 + M|,r|, 



QNc sin2 Ow cos2 Ow Q'^iQ'^ - M|,)2 + M|,r|, 
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(a) 



(b) 



(c) 



Figure 8: qq ^ Z,Z' at LO and NLO (virtual corrections). 






Figure 9: qq Z,Z' at NNLO (virtual corrections). 



,1 Z,l I Z\l Z.l 

9v 9v + 9a 9a 



{Q' - Mi){Q' - M|,) + MzTzMz'Tz' 



96 sin^ Ow cos4 ewNc [{Q^ - M|,)2 + M|,r|,] [(Q2 _ m|)2 + M|r|] 



(29) 



where Nq is the number of colours, T^/^ii is the partial decay width of the gauge boson and the total 
hadronic widths are defined by 



Fz = T{Z hadrons) = ^ r{Z i/jiipi) 

i 

Fz' = T{Z' hadrons) = ^^r(Z' —>■ ipiipi), 



(30) 



where we refer to hadrons not containing bottom and top quarks (i.e. i = u,d,c,s). We also ignore 
electroweak corrections of higher order and we have included the top-quark mass and QCD corrections. We 
have included only tree level decays into SM fermions, with a total decay rate for the Z and Z' which is 
given by 



i=u,d,c,s 

Tz'= Yl ^{Z' ^i^ii^i) + r{Z' ^bb) + T{Z' ^tt) + 3TiZ' -^U) + 3T{Z' ^i^liPl). 



(31) 



Coming to illustrate the contributions included in our analysis, these are shown in some representative 
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(g) (h) (i) 

Figure 10: qq Z, Z' with real corrections at NLO (g) and at NNLO (/i), (i). 




Q) (k) (m) 

(I) 



Figure 11: Anomalous contributions for Z' production in the gg, qq and qg sectors at lower orders. 

graphs. The complete NNLO expressions of the hard scatterings and the corresponding Feynman diagrams 
can be found in [38j. 

• SM QCD contributions 

We show in Fig. [8] the leading 0{aw) and some typical next-to-leading order O(a^as) (LO, NLO) 
contributions to the process in the annihilation channel (virtual corrections). Examples of higher order 
virtual corrections included in the hard scatterings are shown in Fig. O which are of O(a^a^), while the 
corresponding real emissions, integrated over the final state gluons, are shown in Fig. [TU]at NLO (graph g)) 
and NNLO (graphs h) and i)). 

• Anomalous corrections 

We shown in Fig. [11] the leading anomalous corrections to leptoproduction. At ©(a^a^) there is a first 
contribution coming from the interference between graph j) and the leading order qq annihilation vertex 
(graph a) of Fig. [8l The square of the same graph appears in the anomalous corrections at ©(a^a^^). Other 
contributions that we have included are those due to the exchange of a physical axion and goldstone modes, 
which can be removed in the unitary gauge [11] . Of higher order are the contributions shown in diagram k) , 
1) and n), which contribute via their interference with NLO tree level graphs. For instance k) interferes with 
diagram g) of (|10p . while m) interferes with the LO annihilation graph. The analogous contributions in the 
WZ and GS cases are obtained by replacing the triangle graph with the GS vertex, as in Fig. [S] or, for the 
WZ case, with Fig. [3] Notice that in Fig. [S] in the WZ case the anomaly pole is automatically cancelled by 
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Figure 12: Anomalous contributions for the gg gg process mediated by an anomalous Z' at higher perturbative 
orders. 

the Ward identity on the lepton pair of the final state, if the two leptons are taken to be massless at high 
energy, as is the case. Then, the only new contributions from the anomaly vertex that survive are those 
related to the transverse component of this vertex. This is an example, as we have discussed in |12] . of a 
"harmless" anomaly vertex. A similar situation occurs whenever there is no coupling of the longitudinal 
component of the anomaly to the (transverse) external leptonic current. This property continues to hold 
also away from the chiral limit, since the corrections due to the fermion mass in the anomaly have the typical 
structure 

A^upiq, k)_,, = aife^QW^^^ Q - 2m)C^ e^p^^q^k^ + A*™"^ , (32) 

where A*'^""'' is the truly transversal component away from the chiral limit. The most general expression of 
the coefficient Cq is given in Eq.(A.8) of ref. [Uj. Cq is the scalar 3-point function with a fermion of mass 
rrif circulating in the loop. In both mechanisms anomaly (strictly massless) effects are comparable with the 
corresponding contributions coming from the SM for massive fermions. It should be clear by now that in 
the WZ case the anomaly pole is not cancelled, rather an additional exchange is necessary to re-establish 
the gauge independence of the S-matrix (the axion). In DY this sector does not play a significant role due 
to the small mass of the lepton pair. As we have discussed above, the cancellation of the anomaly is due, in 
this case, to the Ward identity of the leptonic current and there is no axion exchanged in the s-channel. 

4.1 Precision studies on the Z resonance 

The quantification of the corrections due to anomalous abelian gauge structures in DY requires very high 
precision, being these of a rather high order. For this reason we have to identify all the sources of in- 
determinations in QCD which come from the factorization/renormalization scale dependence of the cross 
section, keeping into account the dependence on fip and fiji both in the DGLAP evolution and in the hard 
scatterings. The set-up of our analysis is similar to that used for a study of the NNLO DGLAP evolution 
in previous works |48l 149 j . where the study has covered every source of theoretical error, including the one 
related to the various possible resummations of the DGLAP solution, which is about 2 — 3% in DY and 
would be sufficient to swamp away any measurable deviation due to new physics at the LHC. 
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(a) 



(b) 



Figure 13: GS mechanism: anomalous contribution and counterterm for the qq scattering sector. 



These previous studies have been focused on the DY distributions on the resonance peaks, in particular 
on the peak of the Z, where the accuracy at the LHC is of outmost importance for QCD partonometry. 
The presence of anomalous corrections on the Z peak is due both to the anomalous components of the 
Z in the anomalous models and to the interference between the Z' and the Z, that we have taken into 
account. Notice that in DY the treatment of the anomalous corrections to the Z is drastically simplified 
if we neglect the (small) mass of the lepton pair, as usual. In fact, these are due to trilinear (anomaly) 
vertices which involve the BBB, BYY, BW3W3 and BGG gauge fields - in the interaction basis - all of them 
involving interactions of the Stiickelberg field with the corresponding field-strengths of the gauge fields, such 
as hFs A Fg, hFy A Fy, hF\Y A Fy/ and hFc A Fq, where G denotes the gluon field. The only contribution 
that is relevant for the LHC is then one obtained by projecting the hFc ^Fq vertex on the physical axion x, 
whose mass is, in principle, a free parameter of the anomalous models. It is then clear that the axion channel 
plays a more important role in the production of the top, due to its large mass, than in leptoproduction. We 
will now briefly summarize the results for the new contributions in DY, starting from the non-anomalous 
ones. 

In the qq sector we have two contributions involving triangle fermion loops see Fig. [11] k,m. The one 
depicted in Fig. [TTVn is a two-loop virtual correction with a Z or a Z' boson in the final state, while in 
Fig. Illk we have a real emission of a gluon in the final state which is integrated out. The first contribution 
has been calculated in [501 Ell E21 [53] , 

A^g{x,Q^fil,m^) = 6(1 - x)a^'a§CFTf^ (^)' x 

[0{Q^ - 4m2)Gi(mVQ^) + ^(W - Q2)G2(mVQ^)) (33) 

where Cp and Tj are the color factors, q = u,d,c,s, Q = t,b and m the mass of the heavy flavors, while in 
the massless limit the functions Gi and G2 are given by 

Gi{m = 0) = 3 log (^^^ - 9 + 2C(2) 

6*2(771 = 0) = (34) 
and Q represents the invariant mass of the system. The contribution of Fig. Illk in the massless limit is 
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given by 

A« (x,Q2,/x|,m = 0) = a^'a^CpTf^ (^)' x [j^^ ["2 + 2x{l - log(x))]| , (35) 
while in the qg sector we have the contribution shown in Fig. Illi which is given by 
Agg{x,Q^,fil,m^)=afa^TfU^y x [9{Q^ - Am^)Hi{x,Q^ ,m^) + e{4m^ - Q^)H2{x,Q^ ,m^)] (36) 



2 V TT 

with the massless hmit of Hi{x,Q'^,m?) given by 



Hi{x,Q'^,m = 0) = 2x 









2x 


log 











X \ X 



+ 2(1-- 

X 



1 — 2x log I — 



(37) 



Separating the anomaly-free from the anomalous contributions, the factorization formula for the invariant 
mass distribution in DY is given by 

^ = raziQ^Ml) {WziT,Q'') + Wr'^ir^Q^)} 

Wr'^ir^Q^) = Y,f -f,,,(x,/i?j,M^)A-°'"(-,Q2,/i2,) 

^T"{.x, Q\^?p) = A^^ix, Q^fil, m = 0) + A^^{x, Q^,fil,m = 0) + Agg{x, Q^,fil,m = 0) 

(38) 

that we will be using in our numerical analysis below. 
4.2 Di-lepton production: numerical results 

We have used the MRST-2001 set of Pdf's given in [54j and [55j. We start by showing in Fig. [T^ various 
zooms of the differential cross section on the peak of the Z - for all the models - both at NLO and at NNLO. 
We have kept the factorization and renormalization scales coincident and equal to Q, while the mass of 
the extra Z' has been chosen around 1 TeV. The anomaly-free models, from the SM to the three abelian 
extensions that we have considered (free fermionic [40j and U{1)b-l [M] in Fig. [lH while [/(l)g+u appears 
in Tab. ([7]) ) show that the cross section is more enhanced for the MLSOM, illustrated in Fig. [T3^,c. The 
plots show a sizeable difference (at a 3.5 % level) between the anomalous and all the remaining anomaly- free 
models. A comparison between (a) and (c) indicates, however, that this difference has to be attributed to 
the specific charge assignment of the anomalous model and not to the anomalous partonic sector, which is 
present in (c) but not in (a). The anomalous corrections in DY appear at NNLO and not at NLO, while in 
both figures the difference between the SM and the MLSOM remains almost unchanged. 
Moving from NLO to NNLO the cross section is reduced. Defining the iT-factor 

^^^Ifl^^LZ^ ^ K^,^ (39) 
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Figure 14: Zoom on the Z resonance for anomafous Drell-Yan in the /ip = /i/? = Q at NLO/NNLO for all the models. 



in the case of the MLSOM this factor indicates a reduction of about 4% on the peak and can be attributed 
to the NNLO terms in the DGLAP evolution, rather than to the NNLO corrections to the hard scatterings. 
This point can be explored numerically by the (order) variation 



Act 
Act 



Aa (g) (f) + a (g) A(f) 
Wnnlo — cttvloI 



(40) 



which measures the "error" change in the hadronic cross section a going from NLO to NNLO (Act) in terms 
of the analogous changes in the hard scatterings (Aa) and parton luminosities A(j)). The dominance of the 
first or the second term on the rhs of Eq. (j39p is an indication of the dominance of the hard scatterings or of 
the evolution in moving from lower to higher order. The same differences emerge also from Tab. ([5]) and ([6]). 
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da'^'^°/dQ [pb/GeV] for the MLSOM with Mi = 1 TeV, tan/3 = 40, Candia evol. 
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Table 5: Invariant mass distributions at NLO for the MLSOM and the SM around the peak of the Z. The mass of 
the anomalous extra Z' is taken to be 1 TeV with pip — fif; — Q. 

Differences in the resonance region of this size can be considered marginally relevant for the identification 
of anomalous components in this observables. In fact, in [48] a high precision study of this distributions 
on the same peak (in the SM case) shows that the total theoretical error is reasonably below the 4 % level 
and can decrease at 1.5-2 % level when enough statistics will allow to reduce the experimental errors on the 
Pdf's. It is then obvious that the isolation/identification of a specific model - whether anomalous or not 
- appears to be rather difficult from the measurement of a single observable even with very high statistics, 
such as the Z resonance. 

The evolution of the Pdf's has been performed with Candia ^T] which allows independent variations 
oi jjLp and fipi in the initial state. This analysis is shown in Fig. [T^ where we vary /ii? up to 2(3, while we 
have taken 1/2^f ^ f^R ^ 2/xi?. We observe that by increasing both scales there is an enhancement in the 
result and this is due to the logarithms ln//^///|, and ln(5^//u|,, contained in the hard scatterings. The scale 
variations induce changes of about 4% in the SM case at NNLO and about 3.5% in the MLSOM on the peak 
of the Z. Notice that the variations are not symmetric as we vary the scales and the percentual changes 
refer to the maximum variability. This typical scale dependence is universal for all the studies presented 
so far on the peak of the Z and is a limitation of the parton model prediction. After a large data taking, 
optimal choices for the Pdf's and for /x/j and will allow a considerable reduction of this indetermination. 
In subfigure (|15l b) we repeat the same analysis, for the same cm. energy, this time for Q ~ 1 TeV, on the 
Z' resonance in the MLSOM, for a sizeable coupling of the anomalous gauge boson, gs = 1. Compared to 
the value on the Z peak, the reduction of the cross section is by a factor of 2 x 10^. Also in this interval the 
variation of the differential cross section with the two scales is around 3%. We have added a table (TabjT]) 
in which we show results for the total cross sections for the various models at the Z peak. In the first line 
of each column we show the results for the total cross section in [fb], in the 2nd line the total width Tz', 
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dG'^'^^°ldQ [pb/GeV] for the MLSOM with Mi = 1 TeV, tan/3 = 40, Candia evol. 
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Table 6: Invariant mass distributions at NNLO for the MLSOM and the SM around the peak of the Z. The mass of 
the anomalous extra Z' is taken to be 1 TeV with = = Q. 



a?„f° [fb], = 


= 14 TeV, Ml = 1 TcV, tan ,5 = 40 




MLSOM 


f/(l)s-L 




FreeFerm. 


0.1 


5.982 


3.575 


2.701 


1.274 




0.173 


0.133 


0.177 


0.122 




0.277 


0.445 


0.252 


0.017 


0.36 


106.674 


105.567 


53.410 


42.872 




2.248 


1.733 


2.308 


1.583 




4.937 


13.138 


4.991 


0.586 


0.65 


240.484 


143.455 


108.344 


51.155 




7.396 


5.700 


7.592 


5.205 




11.127 


17.853 


10.124 


0.699 


1 


532.719 


317.328 


239.401 


113.453 




17.810 


13.720 


18.274 


12.530 




24.639 


39.491 


22.370 


1.550 



Table 7: Total cross sections, widths and atot x BR{Z II), where BR{Z II) = Tzf^u/Tz'-, for the MLSOM and 
three anomaly-free extensions of the SM; they are all shown as functions of the coupling constant. 
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Figure 15: Zoom on the Z resonance for anomalous Drell-Yan for varying factorization and renormalization scales at 
NNLO for the SM and the MLSOM. Results are shown for Q ~ 91 GeV (a) and 1 TeV (b) both for = 14 TeV. 

expressed in GeV and in the 3rd line the observable atot x BR{Z II), where BR{Z II) = T ^i^ij/T z' ■ 
These quantities refer to the value of the coupling constant Qz listed in the first column. 

We show in Fig. [TBk.b two plots of the results for the MLSOM of the DY cross section on the peak of the 
Z and of the extra Z', where we vary the scales both in the hard scatterings and in the parton luminosities. 
We have added and subtracted the anomalous sector in order to estimate their size respect to the remaining 
contributions. As we have already pointed out, the scale variability at NNLO is larger than the changes 
induced on the result by the anomalous graphs. The anomalous effects are more visible at large Q (subfig. 
(b)), on the resonance of the extra Z' , and are due to the behaviour of the anomalous components at large-x 
due to a growing Q. 

In Fig. [T7] we show in (a) a plot of the MLSOM for different values of gs and for different values of 
and fip. The first peak (purple line) corresponds to qb = 0.1 the 2nd (blue line) to qb = gy and so on. 
As gB grows the width of each peak gets larger but the peak- value of the cross section decreases. Different 
choices oi gB correspond to slightly different values of the mass of the extra Z' because of the relation 
between the Stiickelberg mass Mi and Mz' given in eq. (jl3p . For a fixed value of the coupling, the effects 
due to the variations of the scales become visible only for gs = ^ and in this case they are around 2-3%. 
In the case gs = I (red line), the uppermost lines correspond to the choice fip = 2Q, = ^I'^P^F and 
[IB. = 2/UFi while the lowermost lines correspond to the choice [ip = fiR = 1/2/Xir and = 2^f- Again, 
we notice that if we increase and iir the cross section grows. 

In Fig. 117b we show the result of a comparison between the MLSOM and the anomaly-free extensions. 
We have also included the /j-r/ hf scale dependence, which appears as a band, and the variations with respect 



^Notice that we have chosen = gs for the MLSOM. 
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Figure 16: Plot of the DY invariant mass distributions on the peak of the Z (a) and of the Z' (b). Shown are the 
total contributions of the MLSOM and those in which the anomalous terms have been removed. The variation of the 
result on and is included both in the hard scatterings and in the luminosities {^{^p/ (ir)). 



to gB- As shown in this figure, the red lines correspond to the MLSOM, the blue lines to the U{1)b-l model, 
the green lines to the free fermionic model and the purple lines to C/(l)g+u- Right as before, the first peak 
corresponds to qb = 0.1, the 2nd to gs = gy etc. The peak- value of the anomalous model is the largest of 
all, with a cross section which is around 0.022 [pb/GeV], the free fermionic appears to be the smallest with 
a value around 0.006 [pb/GeV]. 



5 Direct Photons with GS and WZ interactions 

The analysis oipp 77 proceeds similarly to the DY case, with a numerical investigation of the background 
and of the anomalous signal at parton level. 

We start classifying the strong/weak interference effects that control the various sectors of the process 
and then identify the leading contributions due to the presence of anomaly diagrams. 

Direct photons are one of the possible channels to detect anomalous gauge interactions, although, as we 
are going to see, also in this case the anomalous signal remains rather small. Direct photons are produced by 
partonic interactions rather than as a result of the electromagnetic decay of hadronic states. At leading order 
(LO) they carry the transverse momentum of the hard scatterers, offering a direct probe of the underlying 
quark-gluon dynamics. The two main channels in pp collisions are the annihilation qq and Compton {qg), 
the second one being roughly 80% of the entire signal at large pT (from = 4 GeV on). The annihilation 
channel is subleading, due to the small antiquark densities in the proton. The cross section is also strongly 
suppressed (by a factor of approximately 10^^) compared to the jet cross section. For this reason the 
electromagnetic decay of the produced hadrons is a significant source of background, coming mostly from 
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Figure 17: (a) Anomalous Z' resonances obtained by varying gs- (b) Comparisons among anomalous DrcU-Yan in 
the MLSOM versus several anomaly-free models. 
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Figure 18: qq sector for the process qq 77 including virtual corrections at LO (a) and NLO (b,c) 



the vr'^ 77 decay. In this case the angular opening of the two photons is in general rather small, due 
to the small pion mass, and proportional to their energy asymmetry {{Ei — E2)/{Ei + E2)). A limitation 
in the granularity of the detector, therefore, may cause the two photons to be unresolved experimentally, 
giving a spurious signal for direct photons. Even in the presence of enough granularity in the detector, very 
asymmetric decays can also cause the failure of the experimental apparatus to resolve the low energy photon. 
A second source of background is due to ?? ^ 77 decay, which is about 20% of the pion contribution. These 
two contributions account for almost all the background to direct photons. The overall signal to background, 
though small, is supposed to raise with an increasing pT- This is due to the steepening of the tt and 77 spectra 
-at higher Q^- respect to the pT spectra of the parent jet, so there are smaller fractions of these particles that 
can fragment into photons. The reduction of the background can be performed either using reconstruction of 
the contributions due to the pions by measurements of the invariant mass of the pair in the final state, or by 
imposing isolation cuts. In the isolation procedure one can eliminate events with more than 2 particle in the 
final state, considering that emissions from fragmentation are usually accompanied by a large multiparticle 
background. The selection of appropriate isolation cuts are one of the way to render the anomalous signal 
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Figure 19: Real emissions for qq — > 77 at NLO. 




Figure 20: qg sector for the process qg — > 77. 



more significant, considering that the tagged photon signal, although being of higher order in (NNLO in 
QCD), is characterized by a two-photon-only final state. As we are going to see this signal is non-resonant, 
even in the presence of an s-channel exchange, due to the anomaly. 

We show in Fig. [18] a partial list of the various background contributions to the DP channel in 
collisions. We show the leading order (LO) contribution in diagram (a) with some of the typical virtual 
corrections included in (b) and (c). These involve the qq sector giving a cross section of the form 

Gqq = aem(ci C2as)- (41) 

These corrections are the NLO ones in this channel. The infrared safety of the process is guaranteed at the 
same perturbative order by the real emissions in Fig. [19] with an integrated gluon in the final state, which 
are also of 0(ag„as). 

A second sector is the qg one, which is shown in Fig. [20] also of the same order (O(asaem))- These 
corrections are diagrammatically the NLO ones. In general, the NLO prediction for this process are improved 
by adding a part of the NNLO (or O(ag^a^)) contributions, such as the box contribution (j) of the gg sector 
which is of higher order (0(Q!g^a^)) in a^, the reason being that these contributions have been shown to 
be sizeable and comparable with the genuine NLO ones. All these corrections have been computed long ago 
[34j and implemented independently in a complete Monte Carlo in |321 [57] with a more general inclusion of 
the fragmentation. More recently, other NNLO contributions have been added to the process, such as those 
involving the gg sector through Oi^cP'^^oc'g), 

(^gg = al^idial + ^^2"^), (42) 
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Figure 21: gg sector for the process gg — + 77 with virtual and real radiative corrections. 




(W\j 



(o) 




^1 



Figure 22: Anomalous contributions for gg 77 involving the BIM amplitude and its interference with the box 
graphs. 



shown in graphs (A;), (/), (m). The other sectors have not yet been computed with the same accuracy, for 
instance in the qq and qq channels they involve 2 to 4 emission amphtudes which need to be integrated over 
2 gluons. For instance, graph (m) is a real emission in agg which is needed to cancel the infrared/collinear 
singularities of the virtual ones at the same order. 

The anomalous contributions are shown in Figs. [22l [23l [2H [25] and \26^ One of the most important 
partonic process, in this case, is the BIM amplitude shown in diagram (n) which violates unitarity at very 
high energy. In the WZ case, where the axion b appears in trilinear interactions not as a virtual state (such 
as in (bF A F)) this amplitude has two main properties; 1) it is well defined and finite in the chiral limit 
even for on-shell physical gluons/photons and 2) it is non- resonant and can grow beyond the unitarity limit. 
Working in the chiral limit, its expression is given by the Dolgov-Zakharov limit of the anomaly amplitude 
(Eqs. (j2ip and (j22p . with nif = 0), which appears both in the production mechanism of the extra Z', 
{gg — > Z'), where the Z' in the s-channel is virtual, and in its decay into two photons. We recall that the 
presence of an anomaly in the initial and in the final state cancels the resonance of a given channel 
For simplicity we consider amplitude Fig. [22h . assuming to have photons both in the initial and in the final 
state. The amplitude is given by 



. J A r 111 

AbIM = j-^k £[fi,U,h, k2] j^2_M2 



^As stated in the previous sections, the diagram in Fig. 1251 vanish because of a Ward identity. 



29 




Figure 23: Total amplitude for qg 77. 




Figure 24: Another configuration for tlie total amplitude of the qg 77 process. 



= -^e[^i,v,kuk2] [j^j -^ey,u',k[,k'^] 

= -^4f^,^,h,k2] ^ ^e[n',u',k[,k2]. (43) 

where M denotes, generically, the mass of the anomalous gauge boson in the s-channel. If we multiply this 
amplitude by the external polarizators of the photons, square it and perform the usual averages, one finds 
that it grows quadratically with energy. The additional contributions in the s-channel that accompany this 
amplitude are shown in Fig. [271 The exchange of a massive axion (Fig. I27bl. due to a mismatch between 
the coupling and the parameteric dependence between Fig. [27b and b, does not erase the growth (see the 
discussion in the appendix). This mismatch is at the origin of the unitarity bound for this theory analyzed 
in [TI]. The identification of this scale in the context of QCD is quite subtle, since the lack of unitarity in 
a partonic process implies a violation of unitarity also at hadron level, but at a different scale compared to 
the partonic one, which needs to be determined numerically directly from the total hadronic cross section 
app. Overall, the convolution of a BIM amplitude with the parton distributions will cause a suppression of 
the rising partonic contributions, due to the small gluon density at large Bjorken x. Therefore, the graphs 
do not generate a large anomalous signal in this channel. However, the problem of unitarizing the theory 
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Figure 25: Single diagram with an exchanged Z' boson in the s-channel. 




(U) (V) 

Figure 26: Generic representation of the qg 77 process in the presence of a GS vertex of the AVV type. 

by the inclusion of higher dimensional operators beyond the minimal dimension-5 operator bF A F remains. 

The anomalous terms, beside the (n)(n)* contribution with the exchange of an extra Z/Z' which carry 
an anomalous component, which is O(ag^asa^), include the interference between the s/t/u box diagrams 
of gg — > 77 with the same BIM amplitude (n). In the gg sector the anomalous terms give, generically, an 
expression of the form 

'^gg = aem(ai"s"i« + Ci2aia1), (44) 
with the first contribution coming from (n)(o)* and from the interference with the {gg — > 77) box diagram, 
while the second from (n)(n)*. Other contributions which appears at 0(ag^asa«;) are those shown in Fig. [23] 
which involve 2 anomaly diagrams (r) and their interference with the NLO real emission diagram of type 
(m). These contributions are phase-space suppressed. If we impose isolation cuts on the amplitude, we can 
limit our analysis, for the anomalous signal, only to 2-to-2 processes. 

5.1 Helicity amplitudes: massless box diagrams and anomalous interferences 

Moving to the computation of the anomalous contributions to g{pi,±) + g{p2,±) — > 7(/ci,±) + 7(/c2,±), 
coming from the 2-to-2 sector we identify the following non-vanishing helicity amplitudes for the diagrams 
shown in Fig. [22} with the usual conventions 

s = {pi +P2f 

t = {pi- kif = -s/2(l - cos 61) 

u = (pi- kif = -s/2{l + cos e), (45) 



31 



AAA/ 






(b) 



(c) 




{VV\j 



(e) 




Figure 27: Complete list of amplitudes included in the type of graphs shown in Fig. [221 They also have the exchange 
of a physical axion and contributions proportional to the mass of the internal fermion. 



where is the angle between pi and ki , and we obtain 
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where again Z indicates generically either a Z or a Z' and the refers to the contributions with the 
exchange of a axi-higgs. In the above formulas we have omitted ah the couphng constants to obtain more 
compact results. The helicity amplitudes for the massless box contribution have been computed in [58] and 
are given by 
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giving a differential cross section 
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The interference terms are listed in Fig. [22] and the interference differential cross section is given by 
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Here M|°f__ = M^°^++ and we have introduced the parameters of the model, exphcitly given in pT] . 
These correspond to the chiral asymmetries [Dbyy) of the anomalous U{\)b in the BYY mixed anomalous 
diagram whose anomaly is a„; rotation matrix elements from the interaction eigenstates to the physical basis 
after electroweak symmetry breaking (O^); rotation matrix elements in the CP-odd sector (O-^), the strong 
coupling constant g^, the anomalous coupling constant gB and the axial- vector coupling of the neutral gauge 
bosons {Z = Z,Z') to a quark flavour q {g^g)- These parameters combine to define the coupling of the 
physical axion to the various gauge fields iggg,g^-y) 



9gg 



^77 



and 
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(53) 



More details can be found in HD [H [201 [37] . 



5.2 Numerical analysis for direct photons 

In our numerical implementation of double prompt photon production we compare the size of the anomalous 
corrections respect to the SM background evaluated by a Monte Carlo [32l [33]. Since the anomalous 
signal is small compared to that of the SM, we have extracted both for the SM case and the anomalous 
case the gg sector and compared them at hadron level by convoluting the partonic contributions with the 
Pdf's (see Fig. [28]) . In this comparison, the SM sector is given by the graphs shown in Fig. [21] plus the 
interference graphs shown in Fig. [22] In the SM case this second set of graphs contributes proportionally 
to the mass of the heavy quarks in the anomaly loop. At high energy the hard scatterings coming from 
this interference are essentially due to the mass of the top quark running inside a BIM amplitude and are, 
therefore, related to heavy quark effects. In the anomalous case the same set of graphs is considered, but 
now the anomaly contributions are explicitly included. The hadronic differential cross section due to the 
anomalous interactions for massless quarks is given by 
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Figure 28: The NNLO gluon luminosity as a function of r = Q^/S for various value of the invariant mass of the 
photon pair (Q) and energy of the hadronic beam (5) at the LHC evolved with Candia. 
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The contributions which are part of this sector due to exchange of a Z or a Z' and a x (see (a), (b) and (f) 
of the BIM set in Fig. [271) ^^^^ those labelled above, while amt refers to the interferences shown in Fig. 
with the inclusion of a Z' and a physical axion (such as Fig. 127b). 
Defining 

r2-K 
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(55) 



the hadronic cross section takes the form of a product of the gluon luminosity and the partonic gg — > 77 
cross section 

da Q 
dQ~IS' 

5.3 The gg sector 
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Coming to the analysis of the gluon fusion sector, the result of this study is shown in Fig. [29] where we 
plot the gluon contribution to the hadronic cross section for both the SM and the MLSOM, having chosen 
Ml = 1 TeV. We have used the MRST99 set of parton distributions to generate the NNLO gluon luminosity 
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Figure 29: Comparison plots for the gluon sector in the SM and in the anomalous model for a resonance of 1 TeV. 
The box-like contributions are not included, while they appear in the interference with the BIM amplitudes. 



with as{Mz) = 0.1175, Q = ^fiR and = 14 TeV. We have chosen tan/3 = 40 and different values of gB- 
The size of the cross section is around 10~^ [pb/GeV] - right on the mass of the resonance - for both models, 
with a difference that grows as we rise the coupling constant for the anomalous U{1) {gs)- We have chosen 
four possible values for gs'- a small parametric value [gs = 0.1); equal to the coupling of the hypercharge 
gy at the same scale {gs = gy) or to the 31/(2)^ coupling g2 (gs = 52) or, finally, parameterically sizeable, 
with gs = 1- In the interference graphs used for this comparison between the anomalous signal of the 
MLSOM and the SM (in this second case the BIM amplitudes contribute via the heavy quark mass in the 
loops) we have included, beside the BIM amplitude, the entire set of contributions shown in Fig. [27l with 
the exchange of a Z, a Z' and the axi-Higgs x- We have chosen a light Stiickelberg axion with = 30 
GeV. For Q around the anomalous signal grows quite substantially, as we are going to show next. The 
overall flatness of the result in this region - the width of the interval is just 5 GeV in Fig. [22] - shows 
that the corrections are non-resonant and rather small. They are also overlapping for both models and the 
extraction of additional information concerning the anomalous sector appears to be very difficult. Before 
we get into a more detailed analysis of the various contributions to this sector, we mention that we have 
performed isolation cuts on the cross section of the SM background in DP with a choice of i? = 0.4 for 
the radius of the cone of isolation of the photons. This is defined in terms of an azimuthal angle (f) and 
pseudorapidity rj = In tan 0/2, with a maximal value of transverse energy ETmax — 15 GeV in the cone, 
as implemented in [33]. We show a more detailed investigation of the results for the various contributions 
in the gg sector in Fig. [30b . b. The dotted lines are the results obtained by the Monte Carlo and include 
both 2-to-2 and 2-to-3 contributions (pure QCD) with and without cuts, computed at LO and at NLO. The 
size of these contributions is around 2 x 10^^ pb/GeV in the SM case. We show in the same subfigure the 
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Figure 30: (a): SM contributions for the gluon-gluon channel obtained with the Monte Carlo Gamma2MC. These 
are indicated by dotted lines and include all the interferences and the box graphs. Shown are also the anomalous 
contributions of the MLSOM (no box), (b): as in (a) but we have included the Monte Carlo results for the SM qq 
channel at NLO. 

anomalous corrections in the MLSOM, which vary between 10~^ and 10^^ pb/GeV. Therefore, for ~ 1, 
the anomalous sector of the MLSOM (the square of the box terms here are not included for the MSLOM) is 
suppressed by a factor of 10 respect to the signal from the same sector coming from the SM. In subfig. (b) 
we show the same contributions but we include in the SM also the quark channel (shown separately from 
the gluon channel), which is around 10~^ pb/GeV. Therefore, the quark sector overshadows the anomalous 
corrections by a factor of approximately 10^, which are difficult to extract at this value of the invariant 
mass. 

A comparison between the differential cross section obtained by the Monte Carlo and the anomalous 
contributions is shown in Fig. [3Tk . from which one can see that the anomalous components are down by a 
factor of 10'^ — 10^ respect to the background, depending on the value of the anomalous coupling gs- 

A similar comparison between anomalous signal and background for a 1 TeV extra gauge boson but at 
larger invariant mass (2 TeV) of the di-photon is shown in Fig. 131b . In both cases the ratio between the size 
of the anomalous signal and the background is 10~^, showing the large suppression as in the other regions. 

• Anatomy of the gluon sector 

We show in Fig. [32b and b) two plots that illustrate the size of the various contributions to the gg sector 
in the SM and in the anomalous model. We have separated these contributions into several components in 
the MLSOM and SM cases. 

In the anomalous model we have "pure BIM-like" amplitudes: 

1) the square of the "BIM set" shown in Fig. [271 which contains the s-channel exchanges of the Z, of 
the Z' and of the x- Figs. [32k and b these contributions are indicated as "BIM + x". 
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Figure 31: (a) Plot of the invariant mass distribution for Q ~ ITeV and for a Stiickelberg mass of ITeV sliowing 
a comparison between the anomalous contributions (anomalous massless and massive BIM + interference with the 
box) and the SM contributions. All the sectors at 1-loop and 2-loops have been obtained with the Monte Carlo 
Gamma2MC. (b) As in the previous figure but the invariant mass distribution has been chosen around 2TeV. 



Analogously, in the SM case we have that 

1') the BIM amplitudes contribute avi^ay from the chiral limit due to the exchange of the Z gauge boson 
and with top/bottom quarks running inside each of the two loops. In Figs. [32b and b the contributions in 
the SM are denoted by "SM: massive BIM" , and are just obtained by squaring the single BIM amplitude of 
Fig. [27k . In the SM this component is sizeable around the threshold s = ^m1, with mt being the mass of 
the top quark, which explains the cusp in the figure around this energy value. 

The interferences sets of the MLSOM ("BIM-Box") in which the BIM amplitudes interfere with the box 
graphs: 

2) these are denoted as "Z-Box + Zp-Box + x-Box" . The three box diagrams are those shown in Fig. [22| 
corresponding to graphs (o), (p) and (q) in this figure. 

In the SM we have similar contributions. These are generated by 

2') interfering the same box graphs mentioned above with the graph in Fig. [27b . In the SM case, as 
we have already mentioned, these contributions are due to heavy quarks, having neglected the mass of the 
leptons and of the light quarks. This interference is denoted as "Z-Box". 

Finally, 

3) In both cases we have the squared "Box" contributions, which are just made from the o), p) and q) 
diagrams of Fig. [22l 

From a look at Figs. [32b and b it is quite evident that the contributions obtained by squaring the three 
box diagrams are by far the most important at = 14 TeV. We have tried to cover both the region 
Q < 400 GeV and the region 400 < Q < 1000 GeV. The second largest contributions, in these two plots, are 
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Figure 32: Plots of the various anomalous and box-like components in the SM and MLSOM sector at (a) lower Q 
and (b) higher Q. 



those due to the interference between the BIM amplitudes and the box. These contributions are enhanced 
because of the presence of the box amphtude. Notice that in the SM case this interference is small and 
negative for Q < 400 GeV (see Fig. [32k ) and is not reported. In the second region (400 < Q < 1000 GeV) 
this contributions gets sizeable around the tv\fo-particle cut of the scalar triangle diagram due to the top 
quark in the loop and shows up as a steepening in the line labeled as "Z-box" in Fig. 132b . 

The "pure anomalous" contributions, due to the squared BIM amplitudes are down by a factor of about 
10^ respect to the dominant box contributions. A similar trend shows up also in the SM case. Around 
the 2-particle cut the BIM amplitudes both in the SM and in the MLSOM have a similar behaviour, but 
differ substantially away from this point in the larger-Q and smaller-Q region. These are the two regions 
where the effects of the anomaly are more apparent. For instance, for Q < 200 GeV the steep rise of the 
anomalous contributions of the "BIM + x" line is due to the fact that Q is getting closer to the resonance 
of the axion x- In the SM the same region is characterized by a suppression by the ratio nif'^/Q'^. At larger 
Q values, the growth of the anomalous contributions are due to the bad behaviour of the anomaly, which 
grows quadratically with energy, as we have discussed above. This trend is more visible in Fig. 132b (red line) 
in the region Q > 800 GeV. The possibility of increasing the di-photon (gluon-fusion) signal is related to 
the possibility of imposing suitable kinematical constraints in the analysis of the final state - or phase-space 
cuts - in the experimental analysis. To achieve this goal there are several features of the process that should 
be kept into account. The first is the sharp decrease of the cross section at large invariant mass Q of the 
di-photon final state; the second is its increase with the collision energy of the two colliding protons. 
In the first case we have an increase of the relevant parameter r = Q"^ / S characterizing the gluon density; 
in the second case this is reduced considerably, giving a fast enhancement of the gluon luminosity. 

Since bounds on the coupling of the anomalous models may come both from a combined analysis of 
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Figure 33: Double photon invariant mass cross section at Q = Mz plotted as a function of the energy. Here the SM 
contributions include the massive BIM + Box + interference, while the MLSOM contributions include the massless 
and massive BIM + interference. 



both DY and of the di-photon cross section, we have chosen Q around the Z peak (small Q option) and the 
same luminosity evaluated on the peak of the extra Z' (large Q option), both as a function of the collisional 
energy S. This is shown in Fig. [33j The anomalous signal grows with but the non-unitary growth is 
not apparent at LHC energies after convolution with the parton densities; on the other end, the growths in 
the SM and in the anomalous sector of the MLSOM appear to be, at LHC energies, quite similar. 

If an extra Z' is found in DY, then the analysis of di-photon both on the peak of the Z and on the peak 
of the Z' could be used to set reasonable bounds on the coupling of the underlying model, trying to uncover 
the possible presence of an anomalous signal. A r ~ 10^"^ can be reached at the Z peak for a collision energy 
of 10 TeV, which causes an enhancement of the anomalous cross section by a factor of approximately 10^ 
respect to the same anomalous contribution measured for Q ^ 1 TeV. 

In a final figure we show the shapes of the distributions for the SM and the MLSOM using the Monte 
Carlo with a binning of the cross section in both cases. This is provided in Fig. [34l where the colored lines 
on top of each bin indicate the small anomalous signal in the cross section. 

6 Summary and Conclusions 

Both DY and DP have some special features, being characterized by a clean final state. In DY the identifica- 
tion of a new resonance in the neutral current sector would bring to the immediate conclusion that an extra 
Z' is present in the spectrum, but would give not specific indication concerning its true nature. Current 
experimental bounds constrain the mixing of a possible extra neutral component with the Z gauge boson, 
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Figure 34: Binned cross sections for the SM and the MLSOM. The smaU anomalous corrections are the thick hnes 
at the top of each bin. 



with a mass which should be larger than 800 GeV, rendering the future search of extra neutral interactions, 
at least for DY, quite delicate, being the allowed mass range at the tail of the invariant mass distribution of 
this process. For this reason, the identification of a restricted mass range for the extra gauge boson would 
be crucial for its discovery, but unfortunately, there is no extension of the SM that comes with a definite 
prediction for it. 

In intersecting brane models one encounters a similar indetermination and for this reason in our con- 
struction the Stiickelberg mass has been assumed to be a free parameter. If an extra s-channel resonance 
is found, then the investigation of its specific properties would require much more effort and several years 
of data collection at the LHC, especially for larger mass values (above 1.5 TeV or so) of the extra Z'. A 
sizeable width of the resonance could then allow to study the V-A structure of the coupling, though this 
study is quite complex. In most of the models studied so far the corresponding widths are expected to be 
quite small ( ^ 30 — 40 GeV) even for sizeable couplings {gs ~ 0{1)) [40], probably below or barely close 
to the limits of the effective resolution of the detector. 

Under these conditions, deciding over the true nature of the extra Z', whether anomalous or not, would 
then be far more challenging and would require a parallel study of several independent channels. For this 
reason we have analyzed two processes which are both affected by anomalous contributions and could be 
used for correlated studies of the same interaction. 

We have seen that changes in the factorization/renormalization scales both in the hard scatterings and in 
the evolution of the PDF's can easily overshadow the anomalous corrections, making a NLO/NNLO analysis 
truly necessary. We have concentrated our investigation on an extra Z' of 1 TeV in mass and searched for 
anomalous effects in the invariant mass distributions on the Z peak, at 1 TeV and for large Q values (up to 2 
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TeV's). From our analysis the large suppression of the anomalous signal compared to the QCD background 
is rather evident, although there are ways to improve on our results. In fact, our analysis has been based 
on the characterization of inclusive observables, which are not sensitive to the geometric structure of the 
final state. In principle, the shape of the final state event could be resolved at a finer level of detail, by 
analyzing, for instance, the rapidity correlations between the di-photon and a jet, or other similar less 
inclusive cross sections, as in other cases [33]. Obviously, these types of studies are theoretically challenging 
and require an excellent knowledge of the QCD background at NNLO for these observables. This precision 
could be achievable only by a combination of analytical and numerical methods at such perturbative order, 
implemented in a flexible Monte Carlo program. This and other related issues are left to future studies. 
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7 Appendix 



The decay rates into leptons for the Z and the Z' are universal and are given by 
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where i = u,d,c,s and Z = Z, Z'. 

For the Z' and Z decays into heavy quarks we obtain 
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8 Appendix. A comment on the unitarity breaking in WZ lagrangeans 

In the framework of an Abelian AB-model, see [19] for more details, the presence of an untamed growth 
of the amphtude for a 2-to-2 process can be simphfied as follows. In this simple model A is vector-like 
and B is axial-vector like. We also have a single chiral fermion, with an uncancelled EBB anomaly that 
requires an axion-like interaction 6Fb A Fb of the Stiickelberg (6) with the B gauge field for the restoration 
of gauge invariance. The two contributions to BB — > BB are the BIM amplitude see Fig. I35|36l for the 
process mediated by a B boson and a similar one mediated by a physical axion x- If neglect the Yukawa 
couplings (we take the fermion to be massless) the only contributions involved are those already shown in 
Fig- EZl diagrams (a) and (b), but with different couplings. 

^'bi'm = {9BfA'^-^{-h, -fe2)^^f^ L''' - ^] {gBfA>^y-^'{h, k,) (59) 



where Mb = y^Mf + {2gBv)^ is the mass of the gauge boson in the s-channel after symmetry breaking. The 
exchange of the physical axion gives 

B'"'^'"' = ix (^-^aiCsB^ e[fi,u,h,k2]j-^^^e[fi',u',k[,k'2] (60) 

where the overall factor of 4 in front is a symmetry factor, the coefficient ai = comes from the rotation 
of the b axion over the axi-higgs Xi s-nd the coefficient Cbb has been determined from the condition of gauge 
invariance of the anomalous effective action before symmetry breaking 

Cbb - ^4«n^. (61) 
The anomaly diagrams are longitudinal, taking the DZ form 

A'BTM'ian) = {9Bf^{-k')e[f^,.,k^,k,]j^^^ (^/^' - {gBf^{k'')e[^^',u',k[,k',]. (62) 

It is easy to recognize in this anomalous amplitude the same structure of the amplitude for the axi-Higgs 
exchange in Eq. ([60]) . If we add up these two amplitudes and impose the cancellation of the two amplitudes 
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Figure 36: Decomposition of the Stiickelberg axion in a goldstone boson and a physical axion 
(which is a fine tuning) we obtain the condition 

in terms of the mass of the B and the Stiickelberg mass Mi, which docs not have a real solution. This 
simple example that conveys the issue of the presence of unitarity limit for (local) WZ interactions. 
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